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INTRODUCTION 

The physical mechanism responsible for backward scattering of radio waves 
from the middle atmosphere depend on the microstructure of small-scale 
refractivity fluctuations in the vicinity of the Bragg scale, and the spatial 
distribution (or morphology) of this microstructure within the scattering 
volume. The electromagnetic part of the scattering process now appears to be 
reasonably well understood. Characterization of refractivity microstructure, 
however, is only possible in a statistical sense through simplified models for 
turbulent gradient-mixing of passive scalars. Such models were originally 
applied to early tropospheric radio propagation experiments (see e.g. BOLGIANO, 
1968), and with slight modification form the basis for the MST radar technique 
(ROTTGER, 1984). Radar experiments at multiple wavelengths can provide useful 
information on the refractivity microstructure and its dependence on scales 
associated with turbulence (BOLGIANO, 1963; RASTOGI and MATHEWS, 1984). Such 
experiments are feasible only at select facilities. 

In this paper, we discuss the results of a two wavelength (VHF and UHF) 
mesosphere experiment performed at the Arecibo Observatory on January 5, 1981. 

The 46.8-MHz VHF radar (3.21 m Bragg scale) was operated by the Max-Planck- 
Institut fur Aeroncmie (MPI) to provide spectral measurements of signals 
scattered from refractivity fluctuations due to turbulence (ROTTGER et al., 

1983). Other physical parameters such as radial velocities (V r ) , scattered 
signal power (P ) , and Doppler spread (W ) due to turbulence can be derived 
from signal spectra. The 430-MHz UHF radar (0.36 m Bragg scale) was used for 
D-region electron-density (N) measurements using the incoherent scatter 
technique with a comparable height resolution (MATHEWS, 1984). The radars were 
pointed symmetrically about the vertical with a beam spacing of 5.5 degree in 
the meridional plane. Occurrence of a type 4 solar flare during the experiment 
produced enhanced D-region electron-density gradients. This was a unique 
circumstance that provided the possibility of testing the basic premises of the 
turbulent gradient-mixing hypothesis. 

The behavior of physical parameters derived from the VHF experiment is 
analyzed in the next section. We focus on the evolution of a strong layer of 
turbulence at 71 km immediately after the flare onset. It is argued that the 
turbulent layer existed prior to the onset of flare, but was rendered visible 
through enhanced electron-density gradients established after the onset. 
Incoherent-scatter electron-density measurements described in Section 3 clearly 
show substantial enhancement of D-region ionization after the flare. In 
Section 4, we briefly outline the turbulent gradient-mixing hypothesis. The 
VHF and UHF observation are then used to show that the development of the layer 
at 71 km is almost entirely due to enhanced gradients. 
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VHF OBSERVATIONS DURING THE FLARE 

Details of the MPI t V^F experiment at Arecibo and the salient results for 
il^e December** January IftSl^eriod have been reported by ROTTGER et al. (1983). 

We review the observations ‘for January 5, 1981, and present the results of 
further analyses. 

Figure 1 shows contours of constant received signal power as a function of 
height and time. The ST region from 9 to 27 km shows several thin persistent 
layers. The mesospheric signals are weaker in comparison, but layers in the 
vicinity of 65, 71, 73, and 77 km can be distinctly seen. The layer at 71 km 
appears suddenly after 12:16 AST, coincident with the occurrence of the flare, 
but vestiges of it could be seen even earlier. A weak (1-2 dB) but abrupt 
enhancement of received signal power at this time occurs at all heights, and is 
probably due to enhanced solar radio emission seen through a secondary lobe of 
the antenna radiation pattern (ROTTGER, 1983). Details of flare occurrence are 
shown in Figure 2. 

5 JAN 1981 



Figure 1. Contour plots of the total received signal power (P + P ) at 

2 dB increments as a function of height Z and time in AST fuT - n 4 hr) for 
the MPI VHF experiment at Arecibo. Noise enhancement at 12:16 hr AST and 
subsequent development of the echoing region at 71 km is attributed to 
the occurrence of a type 4 Solar Flare. (ROTTGER et al., 1983). 


Figure 3 shows the linearly scaled spectra of the received signal and 
profiles of the noise power P^, signal power P g , radial velocity V r , and 
Doppler spread W g in the mesosphere just before and just after the flare 
onset. A new layer with a near 10 dB enhancement in P appears to form at 71 
km, but the weak layers at 68 and 73 km do not show any Significant variation 
in the signal power. Enhancement of incoming ionizing radiations (Lyman Alpha 
and hard X-*rays) during the flare cannot produce turbulence, but it can amplify 
the refractivity structure within an already turbulent layer by creating 
steeper gradients in the ambient electron densities. 

The intensity of turbulence within a layer of thickness L is characterized 
by the parameter e denoting the rate at which energy is dissipated per unit 
mass by turbulence, e can be related to overturning by buoyancy of the largest 
eddies (of size L) within the layer. Mean background wind shear produces a 
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Figure 2. Solar X-ray fluxes in the 1-8 A band (upper) and 0,5-4 A band 
(lower) recorded by the SMS GOES satellite for January 4-5, 1981. A 
Type 4 solar flare occurred at 16:13 hr UT (12:13 hr AST) on January 
5. It reached a maximum at 16:16:30 hr UT and lasted for 11 min. 
(Solar-Geophysical Data Comprehensive Reports, No. 443 pt II, July 
1981, NOAA) • 


velocity differential across the layer which is mixed by eddies of successively 
smaller sizes. The velocity profile in the vicinity of the layer is slightly 
deformed to maintain a supply of energy to the layer at the rate e. Details 
of this process have been considered by LINDEN (1979) and R0TTGER (1981). The 
distribution of turbulent velocity fluctuations across the layer is mapped into 
the Doppler spectrum of the scattered signal, and the Doppler spread W is 
linearly related to £ . A small contribution to W g due to wind shear ii 
negligible for radar pointing directions close to vertical (SATO, 1981). 

Figure 4 shows the time evolution of received signal power P and Doppler 
spread W g in the vicinity of the 71-km layer. At 70.8 km, estimates of W 
(without any noise correction) are slightly less than 1.5 m/s for noise afone. 
The statistical variation in this parameter is about 0.2 m/s and remains at 
about the same level for time intervals 11:30-12:15 AST and 12:15-13:00 AST. 

We conclude, therefore, that the average energy dissipation parameter <e > 
remains reasonably constant before, during and after the occurrence of the 
flare. 

ELECTRON-DENSITY MEASUREMENTS 

Profiles of electron densities obtained with the incoherent-scatter 
technique are shown in Figure 5. The time resolution for these measurements is 
typically 10 min and height resolution is typically 0.5 km. The profiles have 
been rescaled from total signal power profiles that are smoothed and filtered 
to remove the effect of external interference. A typical uncertainty in 
electron-density measurements is of the order of 10 electrons cm - ^ . The 
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Figure 3. Received signal spectra in the mesosphere, linearly rescaled 
for each height, are shown on the left for the VHF experiment. 
Derived parameters P n , P g , V r , and VJ g (see text) are shown 
in the right panels. Top two panels cover the period just before 
the flare reached its peak. The two lower panels are after the 
flare peak. Appearance of the layer at 71 km is discussed in the 
text. 


profiles are not obtained at equal time spacing. The actual averaging 
intervals are tabulated in Figure 5. Profiles 1 and 2 correspond to periods 
just before and just after the flare onset. Profiles 3 to 5 are for 
successively later periods. At a height of 71 km, electron density abruptly 
increases from -^80 cm" 3 to VOO cm” 3 during the first 10 min of the flare, 
and then gradually decays back to ^100 cm - * over the next half hour. 
Electron-density gradients dN/dZ in the vicinity of the layer show a local 
maximum, and a similar pattern of enhancement after the flare onset. 

DISCUSSION OF TURBULENT GRADIENT MIXING 

The local gradients dN/dZ in electron density N (any passive scalar in 
general) are mixed by turbulence, producing a spatial power spectrum of 
fluctuations in N. The scattered signal is obtained as a component of this 
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spatial fluctuation spectrum evaluated at the Bragg wave number. It follows 
that the scattered signal power P g would depend on [dN/dZ]2. The energy 
spectrum of velocity fluctuations that mix the scalar gradient has several (two 
or more) cutoff scales associated with it. The smallest of these scales 
depends on the energy dissipation parameter e. The exact wave number 
dependence of the energy spectrum, hence the wavelength dependence of P g is 
decided by the location of Bragg scale in relation to the cutoff scales. 

RASTOGI and BCWHILL (1976a, b) discuss wavelength dependence of P g from 
dimensional considerations. B0LGIAN0 (1968) and others have pointed out that 
potential concentrations of the passive scalars should be considered in mixing 
theories. HOCKING (1983) has derived improved relations for wavelength 
dependence of P in the mesosphere for inertial range turbulence. For 
discussion of experiments at Arecibo we consider only the [dN/dZ]^ term, 
since no information can be obtained on wavelength dependence. 

On the basis of earlier discussion, we assume that the turbulent layer 
at 70.8 km has constant thickness L and constant energy dissipation parameter 
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Figure 5. Time evolution of electron-density profiles on January 5, 
1981, using the incoherent-scatter technique at 430 MHz. The 
averaging time corresponding to each profile is shown in the 
table. For the entire time interval shown, the solar zenith 
angle remained within 41.0 to 41.6 degrees. 


e for about 0.5 hr before and 0.5 hr after the flare onset. We will consider 
the following form for the signal power: 

P g = A F(e. L. f Q ) [dN/dZ] 2 + B 

where A is a constant, F contains the dependence on e, layer thickness L, and 
radar frequency f . B is an additive correction term that may be important 
in more exact formulations using potential quantities. We will assume that 
B=0 • 


Figure 6 shows the time variation of the total signal power at 70.8 km in 
the VHF experiment, and electron-density gradients at 2-km scale derived from 
the UHF experiment. Suppose that at two closely spaced times t and t~ 
signal powers are P and P 2# and gradients are [dN/dZ]- and [dN7dZ] 2 . 2 
If all other parameters are constant and B=0, then the ratio R given by: 

R = [P 2 dB - P x dB] / log 1Q {[dN/dZ] 1 / [dN/dZ] 2 } 2 

should be 10. In applying this test to data of Figure 6, we face the problem 
that the power and gradients have disparate time scales. This problem can be 
circumvented by considering transitions in gradients (labeled 1 to 4) and using 
appropriate values of P g , The noise level is reasonable constant at 40 dB. 

For transition 1, we find that R = 8. Transition 2 has only a small change in 
gradient associated with it, so we combine it with 3. These two transitions 
give R = 8.7. Transition 4 cannot be used, as a decrease in gradient is 
associated with enhancement in power. This suggests that the assumption of 
constant £ and L may have become untenable for intervals exceeding 1 hr. 


109 


50i 


P s+n 45^ 


dB 


40 


ifl. 




12 


cm-Ykr 

0 


12 


70*8 km 


M 


n. 


13 hr AST 







1 

1 



3 

4 

y 







13 hr 


Figure 6. Time variation of the total VHF received signal power 
P s*tn for t * ie 71_km la Y er at 3 min intervals (upper panel), 
and the electron-density gradients obtained with the 
incoherent-scatter method with variable averaging intervals 
as shown (lower panel). Gradients are obtained with a 2-km 
scale. Labels 1—4 in the lower panel indicate transitions 
in gradients that are referred to in the text. 


We note that the value of R found in the two cases are both less than and 
close to 10. This suggests that the assumptions made in our analysis hold 
reasonably well, and that enhanced electron-density gradients associated with 
the flare were the principal reason for the layer seen at 70.8 km. A value of 
R>10 is physically unrealistic as it does not allow dependence on other 
parameters. The deficit from R=10 can be attributed to following plausible 
reasons: (a) variability in e, L, (b) local gradients steeper than average 

gradients, and (c) non-zero additive terms derived from mixing of potential 
scalars. In view of intermittent nature of turbulence, we expect (a) to be 
most significant. 
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